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There Is a constant need to he able to solve for enforced 
motion of structures. Spacecraft need to be qualified for ac- 
celeration inputs. Truck cargoes need to be safeguarded from 
road mishaps. Office buildings need to withstand earthquake 
shocks. Marine machinery needs to be able to withstand hull 
shocks. All of these kinds of enforced motions are being grouped 
together under the heading of seismic inputs. 

Attempts have been made to cope with this problem over 
the years and they usually have ended up with some limiting or 
compromise conditions. The crudest approach was to limit the 
problem to acceleration occurring only at a base of a structure, 
constrained to be rigid. The analyst would assign arbitrarily 
outsized masses to base points. He would then calculate the 
magnitude of force to apply to the base mass (or masses) in order 
to produce the specified acceleration. He would of necessity 
have to sacrifice the determination of stresses in the vicinity 
of the base, because of the artificial nature of the input 
forces . 


The author followed the lead of John M. Biggs by ^using 
relative coordinates for a rigid base in a 1975 paper , and 


1. "Introduction to Structural Dynamics" by John M. Biggs, 

McGraw Hill 1964, Sec 6.2. T r 

2. "Fidelity in Shaker Simulation Anaiysis with NASTRAN , T. G. 
Butler, January 1975, Orally presented at the First MSC NASTRAN 

Colloquium. 
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again in a 1981 paper 3 . This method of relative coordinates was 
extended and made operational as DMAP ALTER packets to rigid 
formats 9, 10, 11, & 12 under contract N60921-82-C-0128. This 
method was presented at the twelfth NASTRAN Colloquium. 4 An- 
other analyst in the field, Gary L. Fox, develped a method 5 
that computed the forces from enforced motion then applied them 
as a forcing to the remaining unknowns after the knowns were 
partitioned off. The method was translated into DMAP ALTER' s, 
but was never made operational. All of this activity Jelled into 
the current effort. Much thought was invested in working out 
ways to unshakle the analysis of enforced motions from the limi- 
tations that persisted. In the following theoretical development 
the avenue to complete generality is charted. The method is in 

the process of being cocied and will be implemented as four new 
rigid formats. 

THEORY 


Seismic analysis in the displacement method becomes 
especially challenging, because forces are required in NASTRAN to 
provide loading for the dynamic solutions. The attempt here is 
to admit displacement histories as acceptable loadings by con- 
verting them into equivalent force loadings. The development of 
this theory will start with a statement of the general dynamic 
equation based upon all freedoms being present before any con- 
straints or reductions are applied; this is known as the P-set 


3. "Dynamic Structural Responses to Base Acceleration", Thomas 
G. Butler, Proceedings of the Conference on Finite Element Method 
& Technology, March 1981; Paper No. 8. 

4. "Seismic Analysis Capability in NASTRAN", Thomas G. Butler 
and Robert F. Strang; Proceedings of the 12th NASTRAN Colloquium, 
May 7-11, 1984, pp 92 - 131. 

5. “Solution of Enforced Boundary Motion in Direct Transient and 
Harmonic Problems", Gary L. Fox, Proceedings of the Ninth NASTRAN 
Users Colloquium, Oct 22-23, 1080, pp 96 - 105. 
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(set of ail freedoms obtained from all points, grid and extra) in 
NASTRAN. 


j [«p P ]p 2 + [ b pp> + [ k pp] |{ u P (t >} - { p p (t) }’ m 

where lower case p stands for the differential operator d/dt. 
Freedoms which are directly exposed to seismic forcings (acceler- 
ations, velocities, & displacements) will be given the designa- 
tion "C" (standing for contact freedoms) and the complement of 
this set with respect to the P-set will be designated "J". The P- 
set of Equation (1) will be partitioned between J & C to get 
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Points will be allowed to be loaded with both displacement and 
force histories. This will provide for such cases as a space 
craft being tested in a centrifuge with a shaker on board. In 
such a case there will be body forces being applied by the cen- 
trifuge on ail points including contact points, P c (t), and com- 
plement points, Pj(t); and displacement histories being applied 
by the shaker, u c (t). Single point constraints (SPC's) can be 
applied only to J dof's, out multipoint constraints (MPC's) can 
exist between C & J dof's, however the C freedoms must be chosen 
as independent when defining the constraint. Thus the known 
quantities in equation (2) are the forces on the complement set 
Pj, the forces on the contact set P c , and the displacement his- 
tories at the contact set u c , pu c , and p 2 u c - 

Since the set of u c are known, the terms involving them 
can be expanded from equation (2). Take the known terms in the 


upper partition first: 


j [ M cc ] p2 + [ B ccJ p + [ K cc] { u c} * 
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The dimension of each of these 3 terms is force. Designate the 

C 

set of terms in expression (3) as Pq? i.e. the forces from dis- 
placement histories on the contact freedoms. Next the known 
terms in the lower partition expands into: 


"jcje 2 ♦ 


[ b jc] p + [ k jc] { u c} 


(4) 


the 


Designate the set of terms in expression (4) as Pj? i.e. 
forces on those complement freedoms, J, from displacement his- 
tories due to their being coupled to the contact freedoms, C. 

lcc IP 2 {u c } constitutes forces 
of masses at the contact 

2 K) 

the ac- 


Tne first term of expression (3) [M r 
that develop from the accelerations 

surface. Tne first term of expression (4) |Mj^|p' 
tutes forces tnat develop in the " complement " set from 
celerations of interior masses due to their couplings with the 


consti- 


B ccj p { u c} 


contact set. The second term of expression (3) 
constitutes forces from the speeding of dampers that are con- 
nected between members of the contact set. The second term of 
expression (4) jBj C jpju c j- constitutes forces that develop in the 
" complement " set from the speeding of dampers that are connected 
between the interior and the contact set. The third term of 


expression (3) 
deformation of 


K c dju c j- constitutes forces that develop from the 
elastic elements that are connected between mem- 


(4) 


bers of the contact set. The third term of expression 
jKjcl{ u c} constitutes forces that develop in the " complement " set 
from the deformation of elastic elements that are connected 
between the interior and contact set. The portrayal of the 
forces on the interior dof's must be extracted from the J par- 
titioning of the P-set, otherwise an incorrect distribution would 
result from the increased coupling if they were extracted from a 
reduced order such as N-set or A-set. 
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The scheme here is to treat the excitation histories as 
known only for the purpose of computing forces that develop from 
displacements on contact points. Once the forces from displace- 
ment histories are defined they will be added to boundary force 
histories to give an array of excitations expressed entirely of 
forces in spite of the fact that part develop from displacement 
histories. After the forces from displacement histories are 
fully defined, the contact freedoms u c (t) will henceforth be 
treated as unknown . In effect the scheme is to re-solve for 
displacement histories that are already known. This can be 
characterized with the following example. Put simply; if one 
were to Iook at a single dof system dynamic equation 

rap 2 xlt) + bpx(t) +kx(t)=P(t) (5) 

one could compute the value of the external forcing P(t) if all 
three of the displacement histories were known. For the opposite 
case, one could treat P(t) as known in equation (5), and inte- 
grate it to find the acceleration, velocity and displacement at 
any time. The result would be to recover the values that were 
originally known (assuming perfect differentiation and integra- 
tion routines). This is not an unreasonable approach in view of 
the power in today's computers. 


With the displacements on contact points being treated as 
unknowns, the forces in equation (2) can now be augmented with 
the forces from displacement histories as follows: 
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u_ ( t ) would be recovered if P c <t) & Pj(t) were null. 
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This lays the groundwork for implementation. Provision 
must he made for admitting displacement history specifications as 
hulk data; i.e. p 2 |u(t)j-, p|u(t)J, and u(t). Next, the com- 


putation of P~(t) and Pj(t) must he provided for. Different 
parts of a structure can have certain portions involved in a 
given displacment excitation while other portions could he sub- 
ject to distinctly different excitations. Thus a framework is 
needed for the spatial specification of each distinct excitation. 
There can also he spatially distinct time delays associated with 
individual excitations. But a mechanism already exists in NAS- 
TRAN for such specifications: i.e. DAREA for spatial specifi- 
cation of magnitudes, TABLEDi for time varying amplifications, 
and DELAY for spatial specifications of time delays. All of 
tnese can he used with impunity and without confusion with re- 
spect to the normal input of dynamic data by requring unique set 
ID numbers ana by having a seismic assemoler of enforced load- 
ings. A new case control command called SEISLOAD and a new hulk 
data card called SEISLOAD will be put into service. Bulk SEIS- 
LOAD will act much like TLOADi and RLOAD cards in organizing the 
spatial, temporal, and phase aspects of displacement excitations. 
It will incorporate one additional BCD field to specify the type 
of displacement being input; DISP, or VELO, or ACCE. SEISLOAD 
case control command will activate the hulk SEISLOAD card much 
like the DLOAD case contol command that activates the bulk DLOAD 
card. The Input File Processor (IFP) will assemble the seismic 
bulk data into the initial data block called DYNAMICS. Case 
control will direct the data from its SEISLOAD card to read the 
data from the DYNAMICS data block with a new functional module 
SPD (seismic pool distributor) whose function would be similar to 
the DPD (dynamics pool distributor) to prepare SEISLT (seismic 
load table) and SEISRL (seismic response list) similar to the DLT 
& TRL. Now comes the actual work of processing these tables and 
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lists into actual force histories. SEISLT & SEISRL would be 
input to a second new nodule SEISLG (seisnic load generator) that 
would treat each distinct displacement excitation as an individ- 
ual case. That is, SEISLG would form the partitioning vector of 
the P-set between the C & J sets for one distinct loading. It 
would compute the equivalent set of three force loadings and 
ready it for combining with loads from Load generator modules; 
then turn to the next distinct case and build another partitoning 
vector for this succeeding case and proceed as before in 
computing the equivalent set of three loadings. A record should 
probably be kept for purposes of checking and in setting up 
output sets for recovery of proof of re-solving for the input 
specifications. 

There are several situations that must be anticipated. 
First an important premise must be stated. REGARDLESS OF WHAT 
COMPONENTS OF SEISMIC EXCIATION ARE SPECIFIED (p 2 U, pU, OR U) , 
ALL THREE COMPONENTS EXIST AS A ' CONSEQUENCE OF THE EXISTENCE OF 
ANY ONE OF THEM. For example, if a seismic acceleration were 
given as a specif ication for excitation, the associated velocity 
and displacment histories can be derived by integration. All 3 
components of a seismic disturbance can produce excitation in a 
structure provided that the structure contains appropriate ele- 
ments that are coupled to the contact points. Therefore if only 
one or two out of the three components are specified, the analy- 
sis must be equipped to derive the missing component ( s ) . This 
means that seismic specifications must be differentiated and/or 
integrated to complete the description of the excitation. Modules 
will need to be written to perform both integration and differen- 
tiation of these displacement histories. The options would be 
these when all three components are needed: 


66 


GENERALIZED SEISMIC ANALYSIS 


(a) Only DISP is specified on the SEISLOAD card. 

Consequence: Differentiate twice to obtain seismic velocity 

and seismic acceleration. 

(b) Only VELO is specified on the SEISLOAD card. 

Consequence: Differentiate once to get seismic acceleration. 
Integrate once to get seismic displacement. 

(c) Only ACCE is specified on the SEISLOAD card. 

Consequence: Integrate twice to get seismic velocity and 
seismic displacement. 

Once the three components of seismic excitation are fully enunci- 
ated for one case they will be ready for delivery to SEISLG for 
computation of forces. Each such triplet of histories must be 
identified with its associated spatial companion. Some connec- 
tion must be made with Case Control so as to iceep these various 
combinations of load separated for purposes of managing the 
solution and data recovery operations. 


SEISLG must operate similar to TRLG in that it should 
proauce P-set forces, and D-set forces, and S-set forces. It 
will co this for the C-set based on the SEISLOAD data. It will 
also have to determine which of the J-set are loaded and to what 
extent, due to their individual coupling and prepare these addi- 
tional loadings. After the dynamic load generator has done its 
work on normal forcing, the forces due to displacements should be 
added into the three different partitions of load vectors such as 
the Pp vector. 



+ 




where i represents a distinct contact set. 


(7) 
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For each C dof there exists a distinct set of coupling to 
the J dof 's for mass and for elasticity, and for damping. There- 
fore, for each C dof for each C point there will be a distinct C- 
J partitioning vector. For example, if there are 2 C-points and 
if each point were being excited in 2 translational dof's, there 
are 4 possible couplings for mass, 4 possible couplings for 
damping, and 4 possible couplings for stiffness. Thus there 
would be 3 x 4 = 12 distinct J-C vectors, 12 distinct DAREA 
patterns, 12 distinct TL0AD1 combinations, 2 x 2 X 3 = 12 dis- 
tinct TABLED 1 histories, 3 x 4 = 12 DELAY spatial distributions, 
and 1 SEISLOAD assemblage. 

Translated into a specific example, if the two C-points 
were numbered 50 and 60 and the excitations were in axial (x=l) 
and transverse (y=2) directions, there will be 4 distinct ac- 
celeration histories: 50(x) and 50(y) plus 60(x) and 60(y). The 
mass coupling between 50(x) and its J neighbors would probably 
have a different pattern than that of the mass coupling between 
50(y), 60 ( x ) and 60<y) and their respective J neigbors. So the 
DAREA content for the spatial loading from the acceleration 
excitation at 50(x) will have to be derived from the mass coup- 
ling to 50( x) . Fortunately the DELAY content for the spatial 
time lapse of the acceleration history at 50(x) will be the same 
as the DAREA content for 50(x). Similarly, the DAREA & DELAY 
distributions for 50(y), 60(x), and 60(y) will have to be derived 
from the mass couplings between their J neighbors and at the 
respective points 50(y), 60(x), and 60(y). 

This same pattern of reasoning applies to the formation 
of loadings for displacement histories stemming from stiffness 
coupling between the C dof's and their J neighbors. And again 
this same reasoning applies to the formation of loadings for the 
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velocity histories stemming from damping coupling from the C 
dof's and their J neighbors. TLOADl's and SEISLOAD for the 12 
loadings can be described thusly: 
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DISP 0 50 (y) TLOAD1 
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DISP 0 60 < y ) TLOAD1 
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Now all bookkeeping is in the hands of Case Control and 
the loads are all in terms of force, so the dynamic solution can 
proceed as it normally does, including the recovery of data. The 
output should provide bookkeeping for the several C sets that 
were fed to the SPD (Seismic Pool Distributor module) so that a 
separate reporting of tnese dynamic displacements can be as- 
sembled for comparison with the specified seismic histories 
anc/or a differencing should take place to give a measure of the 
effectiveness in re-solving for the specified seismic inputs. 


APPLICATION 


This theory has been implemented in DMAP form for Direct 
Transients. Although the problems were small pilot examples they 
included extra points and DMIG matrices and involved excitations 
from mass coupling, damping coupling and stiffness coupling. The 
theory has been thoroughly certified. The pilot problem, shown 
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in the plot, represents a simple truss bridge on three founda- 
tions with a seismic wave travelling in the positive x direction 
and disrupting these f oundataions . 



CONCLUSION 

Here at last is an automatic method for handling enforced 
motion that is completely general. The method has been shown to 
be operational in a DMAP mode. There is no special burden on the 
analyst except to provide the usual engineering information 
giving the particulars of his problem. The coding will be com- 
pleted by the summer of 1993 and will be available in the 1994 
release of NASTRAN. 



